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In a two-dimensional parabolic quantum dot charged with A'^ electrons, Thomas-Fermi the- 
ory states that the ground-state energy satisfies the following non-trivial relation: Egs/ihtu) ~ 
N^^'^ fgs{N^^*f3), where the coupling constant, /3, is the ratio between Coulomb and oscillator (hui) 
characteristic energies, and fgs is a universal function. We perform extensive Configuration In- 
teraction calculations in order to verify that the exact energies of relatively large quantum dots 
approximately satisfy the above relation. In addition, we show that the number of energy levels 
for intraband and interband (excitonic and biexcitonic) excitations of the dot follows a simple ex- 
ponential dependence on the excitation energy, whose exponent, 1/0, satisfies also an approximate 
scaling relation a la Thomas-Fermi, Q/{huj) ~ N~''g(N^^^f3). We provide an analytic expression 
for fgs, based on two-point Fade approximants, and two-parameter fits for the g functions. 

PACS numbers: 73.21.La, 68.65. Hb, 73.20.Mf 



I. INTRODUCTION 



Thomas-Fermi theorjfllSlMI has proven to be a valu- 
able tool for the qualitative understanding of atoms and 
molecules. In semiconductor quantum dot^^, which are 
a kind of artificial Thomson atoms with many possibili- 
ties for fundamental research and technical applications, 
Thomas- Fermi theory was shown to agree qualitatively 
and even quantitatively with a mor e e llaborated ap- 
proach like Density Functional Theorjl^, being asymp- 
totically exact in the limit of large electron number^. 

From the computational point of view, Thomas- 
Fermi theory with minor corrections is able to repro- 
duce the ground-state energy of electrons in a quadratic 
potential's at the same level of accuracy of other 
semiclassical or semianalytic approaches like large-D 
expansion^iil or two-point Fade approximant^^^l 

In the present paper, we would like to stress on a less 
studied aspect of Thomas-Fermi theory: the highly non- 
trivial scaling relations following from it. We show that 
the number of electrons, N, and the coupling constant, 
/?, enter the ground-state energy in a scaled form. We 
perform extensive Configuration Interaction calculations 
for quantum dots with 20 < A'^ < 90 in order to verify 
this scaling. In addition, on the basis of the numerical 
results, we show that similar scaling relations are valid for 
the number of excited states in intraband and interband 
excitations. In this way, a universal parametrization of 
the density of energy levels in quantum dots is provided. 

We start with the Hamiltonian of a two-dimensional 
parabolic quantum dot charged with N electrons. In os- 
cillator units, the Hamiltonian can be written as: 



H 



(1) 



The only approximations made in writting Eq. ([T]) 
are the effective-mass description of electrons, the inclu- 



sion of an effective low-frequency dielectric constant, e, to 
model the medium, and the description of confinement by 
means of a harmonic-oscillator potential. These approxi- 
mations are very common and well sustainecP. The cou- 
pling constant (3 = Ecoui/i^^) = e^mi/2/(47rea;^/^/i^/^) 
is the ratio of Coulomb and harmonic-oscillator charac- 
teristic energies. 

The fact that the number of electrons may enter the en- 
ergy in a scaled combination with (3 is, however, not triv- 
ial. Let us write the Thomas- Fermi energy functional^ for 
the present problem: 
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where n{r) is the (surface) density at point r, and a is 
a numerical constant. The above functional should be 
extremized under the constraint 
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Now, it is easy to realize that we can scale r and n in 
such a way that the l.h.s. of Eq. (|3| becomes one, and a 
factor A^-^/^ is extracted from the r.h.s. of Eq. ([2|. As a 
result, we get the following relation for the ground-state 
energy in the Thomas-Fermi approximation: 



(4) 



Notice that the scaled Thomas-Fermi equations depend 
on a single parameter, z = N^^'^f], which combines in a 
particular way the coupling constant and the number of 
electrons. 
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FIG. 1: (Color online) Different contributions to the ground- 
state wave function entering the Configuration Interaction 
calculation. 




II. SCALING IN THE GROUND-STATE 
ENERGY 



FIG. 2: (Color online) Scaling of the ground-state energy 
in medium-sized dots. The large-A'^ expression for the Fade 
estimate, Eq. Kv, is shown as a solid line. 



We first provide an analytical expression for fgg based 
on two-point Pade approximantgl^ in the large- limit. 
It shows that the scaling predicted by Thomas-Fermi the- 
ory is quite general and compatible with true quantum 
effects. 

Let us recall the definition of the ^4,3 Pade approxi- 
mant for the ground-state energy, given in Ref. 1121 which 
interpolates between the — > (perturbation theory) 
and — > 00 (Wigner "crystal") expansions: 



1 + gi/31/3 + 52/32/3 + ■ 



(5) 



We use the large-A^ asymptotic expressions for the 
coefficients^, which lead to the following estimation for 
the ground-state energy: 



Egs ^ 2 



0.698 z+ 1.5 2^/3 + 2.175 z^/^ 
1 + 2.149 zi/3 + 1.5 ^2/3 + 2.175 z ' 



(6) 



In order to verify the universal relation (|6| we per- 
formed extensive Configuration Interaction calculations 
for charged quantum dots. In these calculations, we fol- 
low standard procedures of Quantum ChemistrjElor Nu- 
clear Physics^"'. 

The starting point is the Hartree-Fock solution of the 
problem. Then a basis of functions made up from (i) 
the Hartree-Fock state, \HF), (ii) one-particle one-hole 
(Iplh) excitations, that is \a^) = e\e^\HF) , and (iii) 
two-particle two-hole (2p2h) excitations, i.e. jcrp, /lA) = 
e\e^pe^e\\H F) , is used in order to diagonalize the Hamil- 
tonian. Notice that a < p are single-particle states above 
the Fermi level, and /i < A are states below the Fermi 
level. A schematic representation is given in Fig. [l] In 
the Hilbert subspace with the same quantum numbers of 
the Hartree-Fock state, the electronic Hamiltonian takes 
the form: 



H 




(7) 



where Eh f = 
energy, A^,^,^^^, 
Dankoff matrix. 



HF\H\HF) is the Hartree-Fock 
= {a'pl\H\ap) is the Tamm- 
DHF,app.\ = {HF\H\ap,fj,X), 
Bcr'^L',apti\ = {a'p'\H\ap,pX), and Ccr'p'^'\',appX = 
{a'p',p'X'\H\ap,p,X). D* and B* are, respectively, 
the transposes of matrices D and B. Explicit matrix 
elements are given in Appendix \K\ for completeness. 

In sectors with quantum numbers others than the 
Hartree-Fock state, the first row and column of matrix 
([7]) should be dropped. 

An energy cutoff of 3 hio in the excitation energy is 
used to control the dimension of the Hamiltonian matrix. 
The estimated error in the ground-state energy is below 
0.2 %. 

We computed the ground-state energy of dots with 
N = 20, 30, 42, 56, 72 and 90, and confinement strengths 
huj — 6, 12 and 18 meV. Notice that these are closed shell 
quantum dots with ground-state angular momentum and 
spin quantum numbers L = S — Q. GaAs parameters, 
m — 0.067 Too and e — 12.8, were used. We performed 
the calculations for three-dimensional dots in which the 
confinement along the symmetry axis (the z axis) is mod- 
elled by a rigid-wall well of width, = 25 nm. The con- 
stant NE'f \ where E'f^ = h^n^ /{2mLl), was removed 
from the ground-state energy. The results are depicted in 
Fig. [2] (dots) along with the large- Pade estimate given 
by Eq. (|6| (solid line). Scaling of the ground-state en- 
ergy is apparent. The maximum deviations with respect 
to the Pade estimate are below 10 % for the smallest 
dots with N — 20. Notice that, for the parameters used 
in the calculations, the scaled variable N^/'^f3 takes val- 
ues around 1, i.e. in the transition interval from weak 



3 



to strong coupling^. In order to test the whole interval, 
we use additional control cases: one of them deep in the 
strong coupling regime {N = 42, hio ^ 2 mcV), and the 
other in the weak coupling region {N = 20, hui = 50 
meV). They also fit the Pade estimate. 



III. INTRABAND EXCITATIONS 

We now turn to the intraband excitations. For sim- 
plicity, we consider the excited states of the closed-shell 
quantum dots studied above. We restrict the analysis to 
sectors with the same quantum numbers as the ground 
state, L = S* = 0, in such a way that the ground and 
excited states come out from the same calculation. A 
sample of the results is shown in Fig. [s] (a) for the 42- 
electron dot with confinement huj = 6 meV. First, we 
notice that the excitation gap, which is 2 huj in the non- 
interacting /3 — > limit, is renormalized by Coulomb in- 
teractions to around 6 meV, that is 1 huj. In the oppo- 
site, f3 — *■ oo, limit, the excitation spectrum is that of a 
big (Wigner) molecule, whose small-oscillation frequen- 
cies are independent of The lowest of these frequen- 
cies, i.e. that one determining the gap, should go to zero 
for large N in order to meet the acoustic phonon of the 
Wigner lattice. Then, we can look for a simple interpo- 
lation formula in order to fit the numerical data for the 
excitation gap: 



AEi 



2 + aiP 



(8) 



The parameter 7 appears to be very close to 1/4, thus 
we fixed it to 1 /4 and fit again the data in order to obtain 
fli and bi. The results are shown in Fig. [s] (b) as a 
function of z = N^^^p. We stress that this is only a 
useful representation because AEi does not scales with 
z, even though Fig. |3] (b) shows an approximate scaling 
for intermediate couplings. For the parameters ai and 
61, we get: oi = 3.659, hi = 1.878. The result of the 
fit is excellent, with maximum deviations below 10 %, 
the same as for the ground-state energy. Notice also that 
N^^'^AEi/{Huj) goes to a universal value, ai/61, in the 
strong-coupling limit, /3 — > cxi. Expressions similar to 
Eq. Q for the gap to the first excited state should be 
valid in other angular momentum and spin sectors, and 
also for the energy of collective states (spin- and charge- 
density excitations). 

The second point to notice in Fig. [3] (a) is the ex- 
ponential growth of the number of states for excitation 
energies above 6 meV. This simple exponential depen- 
dence on excitation energy is known in Nuclear Physics 
as the constant temperature approximation (CTA).'i2 



: A^oexp(A£;/e). 



(9) 



^states 

It seems to be a quite general property of the excita- 
tion spectrum of quantum systems. We verified it, for 
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FIG. 3: (Color online) (a) The number of excited states in 
the 42-electron quantum dot as a function of the excitation 
energy. The confinement strength is huj = 6 meV. (b) The 
excitation gap to the first excited state as a function of the 
scaled variable z = N^^^/S. (c) The temperature parameter, 
0, in scaled variables. Fits in (b) and (c) correspond to Eqs. 
@ and |Tol ). 



instance, in the energy spectrum of small quantum dots 
in strong magnetic fieldJI^. 

We fit the numerical data corresponding to the first 25 
excited states of the quantum dots mentioned above in 
order to extract the "temperature" parameter, O, in Eq. 
([9]). We took the first excited state as the reference of 
energy. The next 24 states are only a few (1 to 4) meV 
above the first excited state. 

In order to deduce the universal properties of Q let us 
recall the /? — > and /3 — > 00 asymptotic regimes. In 
the P —f 00 limit, we expect for a behaviour similar to 
AEi, that is N'^^^Q /{hw) should take a universal value. 
On the other hand, in the /3 — > limit the excitation 
energies (with respect to the first excited state) are pro- 
portional to /?, thus we may write a simple interpolation 
formula for the temperature parameter: 
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FIG. 4: (Color online) (a) The interband (excitonic) exci- 
tations in the quantum dot with 42 electrons and hu — 18 
meV. In the x axis the reference energy is the first exci- 
tonic state, (b) Scaling of the temperature parameter for 
the lowest-energy excitonic states. 



Hlu 



b2 + z' 



(10) 



where a2 = 0.360, 62 = 1.226. The quahty of the fit is 
also very good as can be seen in Fig. [s] (c). 



IV. INTERBAND EXCITATIONS: 
STATES 



EXCITONIC 



Next, we study the interband excitonic excitations of 
dots with N = 20, 30 and 42, and fiw = 6, 12 and 18 
meV. The two control cases in the strong and weak cou- 
pling regimes are also included. A basis for excitonic 
states in these dots is build up in the following way: (i) 
states with one additional electron above the Fermi level 
in the conduction band and a hole in the valence band, 
|cr, r) = el^hllHF) , (ii) states with two additional elec- 
trons above the Fermi level in the conduction band, a 
hole in the conduction band, and a hole in the valence 
band, \ap,T,ii) — e\,e^phXe^\HF) . Details of the compu- 
tational scheme in the present case can be found in Ref. 
[T71 The Hartree-Fock single-particle states for holes are 



FIG. 5; (Color online) Structure of the biexcitonic wave func- 
tion in the Configuration Interaction calculations. 



obtained from the Kohn-Luttinger Hamiltonian in the 
presence of the electronic background. In our model cal- 
culations, the oscillator lengths for electrons, heavy holes 
and light holes are equal. Kohn-Luttinger parameters for 
GaAs are usecP^. With a cutoff in the excitation energy 
of 3 hio, the basis dimension is reduced to around 5000. 

We show in Fig. |4](a) a typical spectrum of excitonic 
excitations, corresponding to a dot with = 42 and 
hu) = 18 meV. The states are characterized by the total 
angular momentum J- — + — —3/2, and 

total electronic spin projections, S = 1/2. and Lh are 
orbital angular momenta of electrons and holes, respec- 
tively, and M is the band momentum of holes along the z 
axis. In Fig. |4](a), the x axis excitation energies are mea- 
sured with respect to the first excitonic state. The lowest 
40 states shown in the figure follow two different CTA fits, 
corresponding to AE < 12 meV and 12 < AE < 16 meV. 
The discontinuity in the slope seems to be a quite general 
fact^'' related to different mechanisms of formation of the 
states. 

We use the lowest 10 states in order to find a tem- 
perature parameter in the studied dots, and a law like 
Eq. ( 10 ) to fit the data. The found parameters are: 
03 — 1.373, 63 = 0.559. It is remarkable that the fit per- 
forms very good, as shown in Fig. |4] (b) , signaling that 



the electronic background determines global properties of 
the excitonic excitation spectrum. Unlike the intraband 
excitations, however, we expect the parameters 03 and 63 
to depend weakly on the dot material (GaAs in this case) 
because of the Kohn-Luttinger Hamiltonian entering the 
calculations. We shall test in the future to what extent 
this happens. 



V. INTERBAND EXCITATIONS: 
STATES 



BIEXCITONIC 



Finally, let us consider the interband biexcitonic exci- 
tations in our medium-sized dots. The basis functions 
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FIG. 6: (Color online) (a) The interband (biexcitonic) excita- 
tions in the quantum dot with 42 electrons and hoj = 6 meV. 
In the X axis the reference energy is the first biexcitonic state, 
(b) Scaling of the temperature parameter for the biexcitonic 
excitations. 



for the Configuration Interaction calculations contains 
two additional electrons above the Fermi level in the 
conduction band, and two holes in the valence band, 
\ap,Tri) = el^ejjhlhl^lHF), with a < p and r < ?y. A 
schematic representation is given in Fig. |5] In Appendix 
[B] we give explicit expressions for the Hamiltonian ma- 
trix elements. 

With a cutoff in the excitation energy of 2 hu, the 
Hamiltonian matrix has dimension around 3000. 

We draw in Fig. [6] (a) the spectrum of biexcitonic 
excitations in the dot with 42 electrons and Huj = 6 meV. 
The quantum numbers of the states shown in the figure 
are: = 0, 5 = 0. We see that in a single CTA fit 
we may comprise the first 35 states. These first states 
are to be used in the determination of the temperature 
parameter. 

The scaling of Qxx is shown in the lower panel of Fig. |6] 
Notice the power of N, which is now 1/2 instead of 1/4. 
We verified that, by taking this power, the dispersion of 
points is reduced notably. Thus, we fit the data with the 
formula: 



where 04 = 3.230, 64 = 5.503. The quality of the fit is 
very good. The same comment about the dependence of 
the parameters on the dot material, made above for the 
excitonic states, applies in the present situation. 



VI. CONCLUSIONS 

In conclusion, we have performed extensive Configura- 
tion Interaction calculations for medium-sized quantum 
dots in order to verify universal relations for the ground- 
state energy and the intraband and interband (excitonic 
and biexcitonic) excitation spectrum. The coefficients in 
the r.h.s. of Eqs. ( 6|8|10 ) do not depend even on the 
material the dots are made of. On the other hand, the 
coefficients 03, 63, 04, and 64, we believe, are specific for 
GaAs, but independent of N and huj. 

The work can be extended in many directions. We may 
try to parametrize in a universal way the correlation en- 
ergy, the excitonic and biexcitonic binding energies, the 
excitation gaps to different angular momentum and spin 
sectors, the energy of collective (plasmonic) excitations, 
etc. On the other hand, more efforts towards the un- 
destanding of the empirical relations obtained for the O 
parameters are needed. Research along these lines is in 
progress. 
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APPENDIX A: EXPLICIT MATRIX ELEMENTS 
FOR INTRABAND EXCITATIONS 

In Eq. (|7|, Ehf is the Hartree-Fock total energjiSl 



Ehf — 
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klS^ 



-kls. 
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where pip is the Fermi level, ejf'* is the Hartree-Fock en- 
ergy of the electron state fi, £^°^ is the energy of 2D 
oscillator states, characterized by the quantum numbers 
k (radial number), I (angular momentum), and Sz (spin 
projection). That is: 



&4 + -z ' 



(11) 



41 =^i'^+M2fc +1^1 + 1). (A2) 
The state p is expanded in oscillator states as follows: 
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(A3) 



In the studied closed-shell dots, I and are good quan- 
tum numbers of and the above sum run only over 
k. 

On the other hand, in Eq. ([t]) A is the Tamm-DankofF 
matrixji^ 



+ P{(^',^^\'^/ree\^^',cr), 



(A4) 



where the antisymmetrized Coulomb matrix elements are 
defined as: 



(cr',^|l/ree|^',0-) = (cr', ^| 1/ree ct) - (cr', l/'"ee A^') ■ 

(A5) 

Coulomb matrix elements (ct', /^|l/ree|/^', cr) are com- 
puted in terms of matrix elements among oscillator states 
by using the expansion ( A3 1 . 



Finally, matrices D, B and C are explicitly written as: 

DHF.appX = '^I l/''ee IpTo^) • (A6) 

\\l/ree\^^' ,p)5 



B. 



a' fj.' ,(Tpfj,X 



' ee 

{(j',X\l/ree\p^)Spp' 
{a',ll\l/ree\^)6xp'} . 



(A7) 



C a' p' ^' \' .crp^X 



aa' ^ pp' 



+ (p' , A 1 1 /ree I A' , p) (5o-ct' + (p' , ^ 1 1 /ree I p. A') ^crcr' <5ap' 

+ {p' ,\\l/ree\p,p.')5aa'5p\' + {p' , p\l /ree\p' , p)5aa' 5xx' 
-|-(p', A|l/ree|CT, X')5pp'5p„- + (cr', A|l/ree|p, A')(5^^'(5o-p/ 
-|-((t', \\l/ree\X\<y)5pp'5pp, + (p',Cr'|l/ree|pro^)<^^^'<^AA' 
+ (p',A'|l/''ee|A',Cr)JAM'4'T' + {(^\p\^/ree\^' ,P)5\p'5„p> 

+ {a',p\l/ree\a; \')Sxp'Spp' + (p', A|l/ree|p', cr)(5^A'^p<T' 

-f (cr', A|l/r ee\<^T P')Sfj,X'Spp' + (cr , A|l/ree|p', p)(5^,A"^(Tp' 
+ (p',A^|l/ ''ee|c, p')(5AA"5pcr' + (cr , /ijl/reelp, p')(5AA"^crp' 

-f'((T',p|l/ree|pV7)(5AA"5pp4- (^8) 



APPENDIX B: EXPLICIT MATRIX ELEMENTS 
FOR BIEXCITONIC EXCITATIONS 



In the biexcitonic sector, the Hamiltonian matrix ele- 
ments take the form: 



{a'p',T'r,'\H\ap,Tr^) = 

(EhF + 4'^ + 4"^ + 4"' + 4''^) 5.a'5pp,5rr'5^ 
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+ {a,T'\l/reh\Q.,T)Srjrj']Saa'Sp 



{a,r]'\l/reh\a,T)Sr, 



(Bl) 



The total Hamiltonian, H, in addition to the terms of 
Eq. ([ij, now includes the single-particle energy of holes, 
electron-hole, and hole-hole interactions. The Hartree- 
Fock electron and hole states should be expanded in os- 
cillator functions when Coulomb matrix elements are to 
be computed. For hole states, coming from a Kohn- 
Luttinger Hartree-Fock problem, we have the expansion: 



E 



(B2) 



where k and I are oscillator quantum numbers, m = 
±3/2, ±1/2 is hole (band) angular momentum projec- 
tion, and fcz = 1, . . . , 6 labels sub-band states in the well. 
The relatively large number of terms entering the ex- 
pansion ( B2 1 makes the calculation of Coulomb matrix 



elements for holes lengthy. 
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